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ABSTRACT 

^ : 
o , 

i We investigate the relationship between the mass and velocity fields of the intergalac- 

tic medium (IGM) and dark matter. Although the evolution of the IGM is dynamically 
Ci ' governed by the gravity of the underlying dark matter field, some statistical properties 

^ ■ of the IGM inevitably decouple from those of the dark matter once the nonlinearity of 

ly^ ' the dynamical equations and the stochastic nature of the field is considered. With sim- 

. ulation samples produced by a hybrid cosmological hydrodynamic/N-body code, which 

y—i . is effective in capturing shocks and complicated structures with high precision, we find 

i that the one-point distributions of the IGM field are systematically different from that 

I of dark matter as follows: 1.) the one-point distribution of the IGM peculiar veloc- 

I ity field is exponential at least at redshifts less than 2, while the dark matter velocity 

' field is close to a Gaussian field; 2.) although the one-point distributions of the IGM 

■ and dark matter are similar, the point-by-point correlation between the IGM and dark 

■ matter density fields significantly differs on all scales and redshifts analyzed; 3.) the 
Q^i one-point density distributions of the difference between IGM and dark matter fields are 
Q ! highly non-Gaussian and long tailed. These discrepancies violate the similarity between 
,H I the IGM and dark matter and cannot be explained simply as Jeans smoothing of the 

IGM. However, these statistical discrepancies are consistent with the fluids described 
by stochastic-force driven nonlinear dynamics. 
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^ • Subject headings: cosmology: theory - large-scale structure of the universe 



1. Introduction 

The mass field of the universe is dominated by dark matter with only a tiny fraction of cosmic 
matter in the form of baryonic particles. Gravitational clustering of dark matter is the major 
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imputes for cosmic structure formation. Collapsed dark matter halos host various baryonic light- 
emitting and absorbing objects and information about dark matter is available only via baryons. 
The major component of baryonic matter is in the form of gas, the intergalactic medium (IGM) 
(Here we use IGM to mean baryonic gas both before and after the galaxy formation). Therefore, 
the dynamical relationship between baryonic gas and the underlying dark matter field is essential 
in order to understand the origin and evolution of cosmic structures. 

Since the IGM is only a tiny fraction of cosmic matter, it is usually assumed that its clustering 
behavior on scales larger than the Jeans length traces the underlying dark matter field. The 
density and velocity distribution of the IGM is considered to be the same as the dark matter field 
point-by-point on scales larger than the Jeans length scales. That is 

Sigm{y^, t) = (^dm(x, t), Vig„i(x, t) = Vd^(x, t) (1) 

where 6igm and Sdm are the mass density contrasts and Vigm and v^^ are the velocity fields, 
smoothed on the Jeans length scales. Eq.(l) applies during the linear regime. Even if the IGM 

is initially distributed differently than the dark matter, linear growth modes will lead to eq.(l) on 
scales larger than the Jeans length (Bi, Borner & Chu 1992, Fang ct al. 1993; Nusser 2000, Nusser 
&: Hachnclt 1999, sec also Appendix B). In this case, all statistical properties of the IGM at scales 
greater than the Jean's length are completely determined by the dark matter field. In other words, 
the dynamical behavior of the IGM field can be obtained from the dark matter field via a similarity 
mapping (e.g., Kaiser 1986). 

Because the nonlinear evolution of the IGM is also driven by the gravity of the underlying dark 
matter, it is often assumed that the relation (1) will hold even in the nonlinear regime. However, 
this assumption is probably not valid. It has been shown in hydrodynamic studies that a passive 
substance generally decouples from the underlying field during nonlinear evolution (for a review, 
Shraiman &: Siggia 2001). For instance, a passive substance might be highly non-Gaussian even 
when the underlying field is Gaussian (Kraichnan 1994). This non-linear decoupling is generic to 
systems consisting of a "passive substance" and an underlying stochastic mass field. 

The IGM and dark matter fields interact stochastically since the initial perturbations of the 
cosmic mass and velocity fields are randomly distributed. Their mass and velocity fields act as 
random variables. The importance of the stochastic nature for the evolution of cosmic mass and 
velocity field has been emphasized in many studies (e.g. Bcrcra & Fang 1994, Jones, 1999. Buchert, 
Dominguez &: Peres- Mercader 1999, Coles & Spencer 2003, Ma & Bcrtschinger 2003). In this paper, 
we will address the statistical discrepancy between the random fields of the IGM and dark matter 
in nonlinear regime. 

Statistical discrepancies between the IGM and dark matter mass fields may have already been 
detected in previous studies. The adhesion model of the IGM shows that the deterministic relation 
5b(x, t) = (5(x, t) no longer holds (Jones, 1999). The stochastic nature of cosmic mass field indicates 
that the IGM should be described by a random-force driven Burgers' equation (Matarrese & Mo- 
hayee 2002) which does not always yield the deterministic solution 5b(x, t) = (5(x, t). Cosmological 
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hydrodynamic simulations also find that pigm does not tightly correlate with p^^, but is largely 
scattered around the line 5igm = Sdm, with this scatter not due to noise (Gnedin &; Hui 1998). It 
has also been found that the scatter defined by 

A(5(x, t ) = Sdm{^,t)- Sigm (x, t) (2) 
Av(x, t) = Vdmi^, t) - Yigmi'X-, t) 

is highly non-Gaussian (Feng, Pando &l Fang 2003). This strengthens the conclusion that the 
discrepancy between 5igm and ddm is not due to computational noise or other Gaussian processes, 
but probably arises from the nonlinear evolution of the random fields. 

Here we study this discrepancy at a more fundamental level. As a first step, we concentrate 
on the density and velocity field one-point IGM distributions, and their discrepancy from dark 
matter. The outline of this paper is as follows. §2 addresses the dynamical mechanisms, and gives 

predictions on the statistical discrepancy between the IGM and dark matter density and velocity 
fields. §3 presents the cosmological hydro simulation scheme and samples. The relevant one-point 
statistics are developed in §4. §5 shows the discrepancy and tests the predictions with one-point 
statistics of the density and velocity fields. Finally, conclusions and discussions are in §6. 

2. Statistical discrepancy between the IGM and dark matter fields 

During the linear regime, thermal diffusion leads to a discrepancy between the IGM field and 
dark matter field on scales up to the Jeans length. That is, the IGM mass density perturbations 
are suppressed on scales less than the Jeans length. The IGM density perturbations on scales 
larger than the Jeans length are Gaussian if the dark matter density perturbations are Gaussian. 

However, in the nonlinear regime, a statistical discrepancy very different from a simple thermal 
diffusion can appear. This can be illustrated by considering the isothermal model of the IGM. In 
this case, the IGM density is given by Pigm{^) oc exp[—m(j){x)/kBT], where ^(x) is gravitational 
potential, T local temperature, m the mass of the IGM particles and the Boltzmann constant. 
If (/)(x) is a Gaussian random field, pigm will be a lognormal random field (Zeldovich, Ruzmaikin & 
Sokoloff 1990) and a statistical discrepancy between pigmi^) and dark matter arises. Although the 
isothermal model is not realistic on large scales in general, it reveals that the statistical discrepancy 
between the IGM and dark matter might arise if 1.) the IGM evolution is nonlinear and 2.) the 
relevant fields, like 0(x), are stochastic. 

2.1. The statistical discrepancy of velocity fields 

A more realistic discrepancy appears when considering the differences between the dynamics 
of the peculiar velocity fields of the IGM Vigm and dark matter \dm- Since dark matter particles 
are collision-less, the intersection of the dark matter particle trajectories will lead to a multi-valued 
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velocity field. On the other hand, as a fluid, the velocity field of the IGM will always be single- 
valued. At the intersection of dark matter particle trajectories, the IGM velocity field will be 
discontinuous and yield shocks or complicated structures (Shandarin &: Zeldovich 1989). Shocks in 
the IGM can significantly change the mass density and velocity of the baryon gas, but will exert 
no direct effect on the dark matter field. It is at this point that the dynamical similarity between 
the IGM and dark matter is broken. 

We analyze this situation using the dynamical equations for dark matter and the IGM. For 
growth modes, the peculiar velocity field of the dark matter is vortex-free. One can define a velocity 
potential if as Wdm = —{^/o)^Vdm-, where a is the cosmic scale factor. The dynamical equation of 
the velocity potential is (see Appendix §A) 

^__L,v,,„)-,, (3) 

where (j) is the gravitational potential and depends on density perturbation 5drn = [Pdrn ^'PdnM'Pdrn 
via the Poisson equation (A3). Generally, the field is Gaussian, or only slightly deviates from a 
Gaussian field with (0) = and variance {(jP'). Equation (3) is valid till the intersection or shell 
crossing of dark matter particle trajectories has occurred. 

For the IGM growth modes, one can also define a velocity potential i^igm , by Vjg^ = — ( 1 / a)Vipigm ■ 
The dynamical equation (3) are approximately given by (see Appendix §B) 

^ - ^,i^^^9mf - ^VV..^ = <P. (4) 

^ lxmp{d\B.D{t)/dty ^ ' 

where D{t) describes the linear growth behavior. The diffusion term v in eq.(4) is given by the 
Jeans smoothing. The wavenumber of the comoving Jeans scale is k'j = {a^ /t^){fimp/^kBTo), 
where nip is proton mass. The parameters 7 and /x are, respectively, the polytropic index and 
molecular weight of the IGM. 

The nonlinear equation (4) is actually the stochastic-force driven Burgers' equation or the 
KPZ equation (Kardar, Parisi Sz Zhang 1986; Berera &; Fang 1994). Fields governed by eq.(5) have 
been extensively studied to model structure formation (Barabasi Sz Stanley 1995). The solution to 
eq. (4) depends on two characteristic scales: 1.) the correlation length Tc of the random field (p; 
2.) the dissipation length d = 

^3/4^1/2 

(0^) ^/^{d\nD/dt) which is due to Jeans smoothing. 
The behavior of the field ipigm governed by eq.(4) is determined by its Reynolds number defined as 
7^ = {vc/d)^/^ (e.g., Lassig 2000) or 

^ = (W/'(^)'^'(<^LW)^/' (6) 

where Vc is the comoving correlation length, and Sdmik) is the Fourier component of the density 
contrast on wavenumber k. To derive eq.(6), we assume that the gravitational potential <p is only 
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given by the dark matter mass perturbation. When the Reynolds number is larger than 1, Burgers 
turbulence occurs in the figm field. 

The correlation length Tc of the gravitational potential ^ is larger than the Jeans length and we 

have kjrc > 1- Therefore, TZ can be larger than 1 on scales larger than the Jeans length, even when 
^dmikj) is on order 1. That is. Burgers turbulence develops in the IGM ifigm field while the dark 
matter mass density perturbations are still quasi-linear or weakly nonlinear. We can conclude that 
the evolution of the PDFs of Vigm{^,t) and v<i^(x, t) should be significantly different, the former 
becoming non-Gaussian earlier than the latter. 

When Burgers turbulence develops in the IGM, (figm is characterized by strong intermittency 
and contains discontinuities, or shocks. The probability distribution function (PDF) of ipigm is long 
tailed (Lassig 2000). The intermittent spikes are the events that populate the long tail. This feature 
is supported by observations of the Lya flux transmission (Jamkhedkar, Zhan & Fang 2000, Pando 
et al. 2002, Jamkhedkar ct al. 2003). That is, the randomly distributed shocks and intermittent 
spikes lead to the statistical discrepancy between the velocity fields of the IGM and dark matter. 
One can then expect that the PDF of the difference, Av(x), is long tailed. Moreover, the non- 
trivial difference between the IGM and dark matter holds on the scales larger than the Jeans length 
of the IGM. 



2.2. The 2"°' moment of the density distributions 

The non-trivial velocity difference Av(x, t) will lead to a non-trivial density difference between 
the IGM and dark matter. The linearized continuity equations for dark matter [eq.(Al)] and IGM 
[eq.(Bl)] are, respectively 



dt a 



+ -V-Vd„ = (7) 



^ + iv.v,,„ = 0. (8) 
Thus, using the definition for A(5(x, t), we have 

A(5(x) = 5dm - ^igm = - j dtV ■ Av. (9) 

This result shows that the IGM mass field 5ig„i{^) does not trace the dark matter field ddm{^) 
point-by-point due to the discrepancy in their velocity fields. 

The discrepancy in the mass fields can be measured by the 2"'^ moments of the density fields 
defined as 

_([A^(x)nv2 

_ ([A(5(x)]^)V^ 
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and 

^ ([A^(x)]^)V2 

If (^iprn(x) is perfect tracer of 5dm (x), we have Idm = hgm = 1 = 0. If both (^igm(x) and are 
statistically independent, we have 1 = 1. Thus, from eq.(9), we expect that /, Idm o-'^d ligm are 
smaller at higher redshift, and larger at lower redshifts. 



2.3. The density PDF discrepancy 

Prom the continuity equations eqs.(Al), eq.(Bl) and the definitions eq.(2), we have 

= _iv • [Av + iAd)vam + iAv)digm] . (13) 
We define a window sampling given by 

A«(xo) = J Wr{^' - xo) A 5(x')dx', (14) 

where the normalized window function Wji{'x' — xq) is non-zero around xq on spatial scale R. 
Since AS and ASigm are statistically isotropic, for an isotropic window function WOj(x') the terms 
containing V A (5 and VSigm are negligible. Thus, eq.(13) becomes a Langevin equation for A^(x) 

^ = -iAA«A. (15) 

at a a 

The first term on the r.h.s. of eq.(15) is a friction term with friction coefficient given by 

X = gw J Wnix' - xo) (Vx' • Vdm)dx', (16) 

where the factor gw depends only on the window function (see §5.3). The random driving force rj 
of eq.(15) is 

r) = -J d^'Wni^' - xo) (V • Av)(l + digm). (17) 

Since \Sigm\ < 1 at most places, we can drop Sigm in eq.(17). 

Eq.(15) contains an additive stochastic force, 77, and a multiplicative stochastic force, A. Both 
stochastic forces are given by the random velocity field. The Langevin equation eq.(15) with both 
additive and multiplicative noise is typically used to model intermittent fields (Graham, Hohnerbach 
& Schenzle, 1982; Piatt, Hammel & Heagy 1994, Nakao 1998). In the nonlinear regime, most 
locations in the dark matter density field have ddm < 0, and therefore, V • \'dm and A are larger 
than zero. The friction term A leads to the decay of A^. That is, the density perturbation difference 
A^ caused by the noise 77 tends to zero on average. Although A is positive in most cases, it can 
go negative due to fluctuations. In these cases, A^ will be amplifled exponentially and attain 
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large values, which are the spikes in the field. This is intermittency and the PDF of A will be 
long-tailed. It has been shown that when the stochastic terms ij and A are Gaussian, the PDF 
of is a power law in general (Appendix C). Although 7] and A for the dark matter given by 
eqs.(16) and (17) are not Gaussian, we still can conclude that the PDF of A^ is generally long 
tailed, as long tails are a common result of a multiplicative stochastic force. For instance, if one 
ignores the additive noise term rj, the solution for eq.(15) is A^ oc exp(— / ^\dt). Thus, the PDF 
tail A''^ is longer than that of A so that, for instance, when A is Gaussian, A^ is lognormal. Using 
Appendix eq.(C7), wc sec that the PDF of A^ is generally highly non-Gaussian on scales larger 
than the Jeans length. It is flat in the central part and gradually becomes a power law with index 
not lower than -1. 

By closely studying the dynamics of the velocity fields we have found that the statistical 
discrepancy between the dark matter and IGM velocity fields will manifest itself as follows: 

1. The evolution of the PDF's of Vjg^ and v^m will be significantly different. The former becomes 
non-Gaussian earlier than the latter. Further the PDF of the difference Av(x) will be long 
tailed. 

2. The second moments /, Idm^ and ligm (eqs 10-12) will be smaller at higher redshifts and larger 
at smaller redshifts. 

3. The PDF of A^ are generally highly non-Gaussian on scales larger than the Jeans length. 
In the following sections we will test these predictions 

3. Hydrodynamic simulations 

3.1. WENO hydrodynamic simulations 

To simulate the IGM, we use the the hydrodynamic equations of the IGM in the form of 
conservation laws, eqs.(B4)-(B6). Although the momentum equation (B5) is the typical Navier- 
Stokes equation, gravitational instability leads to a system dominated by growth modes and the 
dynamical equations essentially become a Burgers equation if only the growth modes are considered 
(Berera & Fang 1994). It is well known that the Burgers equation does not reduce initial chaos, 
but increases it (Kraichman 1968). That is, when the Reynolds number is high, an initially random 
field always yields a collection of shocks with a smooth and simple variation of the field between the 
shocks. We conclude that an optimal simulation scheme should capture shock and discontinuity 
transitions as well as to calculate piecewise smooth functions with a high resolution. 

For these reasons we do not use schemes based on smoothed particle hydrodynamic (SPH) 
algorithms. It is well known that one of the main challenges to the SPH scheme is how to handle 
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shocks or discontinuities because SPH schemes smooth the fields. This problem is not yet well 
settled (e.g. Borve, Omang, & Trulsen, 2001, Omang, Borve, & Trulsen 2003). Instead, wc will 
take an Eulerian approach to simulating the IGM. However there is a basic problem in Eulcrian 
based codes in that they cause unphysical oscillations near a discontinuity. An effective method to 
reduce the spurious oscillations is given by designed limiters, such as the total- variation diminishing 
(TVD) schemes (Harten, 1983). However, TVD accuracy degenerates to first order near smooth 
extrema (Godlewski & Raviart 1996). This problem is serious in calculating the difference between 
fluid quantities on the two sides of the shock when the Mach number of a gas is high. But this is 
exactly the case in the gravitationally coupled IGM and dark matter system. In this system, the 
IGM temperature is generally around 10^^^ K and the sound speed a few km s~^ to a few tens km 
s~^, while the IGM rms bulk velocity is on order of hundreds km (Zhan &; Fang 2002). Hence, 
the Mach number of the IGM can be as high as ~100 (Ryu et al 2003). 

To overcome this problem, two algorithms ENO and then WENO were developed (Harten et 
al. 1987, Shu 1998, Fedkiw, Guillermo & Shu 2003, Shu, 2003.) The precision of the WENO scheme 
is found to be higher than that of TVD and the piecewise parabohc method (PPM) for both strong 
and weak shocks. TVD schemes degenerate to first-order accuracy at locations of smooth extrema, 
while the ENO and WENO schemes maintain their high-order accuracy. Both ENO and WENO 
use the idea of adaptive stencils in the reconstruction procedure based on the local smoothness of 
the numerical solution to automatically achieve high order accuracy and non-oscillatory properties 
near discontinuities (Liu, Osher, &; Chan 1994; Jiang Sz Shu 1996). This scheme is probably the 
first successful attempt to obtain a self-similar (no mesh size dependent parameter), uniformly high 
order accurate, yet essentially non-oscillatory, interpolation for piecewise smooth functions. WENO 
is robust and stable and simultaneously provides high order precision for both the smooth part of 
the solution and sharp shock transitions. 

WENO has been successfully applied to hydrodynamic problems containing shocks and com- 
plex structures, such as shock- vortex interaction (Grasso & Pirozzoli, 2000a, 2000b), interacting 
blast waves (Liang &; Chen, 1999; Balsara Sz Shu 2000), Rayleigh- Taylor instability (Shi, Zhang & 
Shu 2003), and magneto-hydrodynamics (Jiang &; Wu, 1999). WENO has also been used to study 
astrophysical hydrodynamics, including stellar atmospheres (Zanaa, Velli &: Londrillo, 1998), high 
Reynolds number compressible flows with supernova (Zhang et al. 2003), and high Mach num- 
ber astrophysical jets (Carrillo et al. 2003). In the context of cosmological applications, WENO 
has proved especially adept at handling the Burgers' equation (Shu 1999). Recently, a hybrid 
hydrodynamic/N-body code based on the WENO scheme was developed and passed typical reli- 
ability tests including the Sedov blast wave and the formation of the Zeldovich pancakes (Feng, 
Shu & Zhang 2004). This code has been successful in producing the QSO Lya transmitted flux, 
including the high resolution sample HS1700+6416 (Feng, Pando & Fang 2003). The statistical 
features of these samples are in good agreement with observed features not only on second order 
measures, like the power spectrum, but also to orders as high as eighth order for the intermittent 
behavior. The code also has been shown to be effective in capturing gravitational shocks during the 
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large scale structure formation (He, Feng and Fang 2004). Hence we believe that for the purposes 
of this work the Eulerian code based on the WENO scheme is the best approach. 



3.2. Samples 

For the present application, we run the hybrid N-body/hydrodynamic code to trace the 
cosmic evolution of the coupled system of dark matter and baryonic gas in a flat low density 
CDM model (ACDM), which is specified by the cosmological parameters {Q,rn,^A, h, as,^b) = 
(0.3,0.7,0.7,0.9,0.026). The baryon fraction is fixed with the constraint from primordial nucle- 
osynthesis as Qi, = 0.0125/i~^ (Walker et al. 1991). The linear power spectrum is taken from the 
fitting formulae presented by Eisenstein &; Hu (1998). 

Atomic processes including ionization, radiative cooling and heating are modeled similarly as 
in Cen (1992) in a plasma of hydrogen and helium of primordial composition {X = 0.76, Y = 0.24). 
Processes such as star formation, and feedback due to SN and AGN activities are not taken into 
account as yet. A uniform UV-background of ionizing photons is assumed to have a power-law 
spectrum of the form J(i^) = J21 x 10~^^(i^/^'//7)~"erg s~^cm~^sr~^Hz~^, where the photo-ionizing 
flux is normalized by the parameter J21 at the Lyman limit frequency i^hi, and is suddenly switched 
on at z > 10 to heat the gas and re-ionize the universe. 

The simulations are performed in a periodic, cubic box of size 25 h~^Mpc with a 192^ grid 
and an equal number of dark matter particles. The simulations start at a redshift ^; = 49 and the 
results are output at redshifts z =4.0, 3.0, 2.0, 1.0, 0.5 and 0.0. At each output stage, we read pdm, 
Pigm (both are, respectively, in the units of p^^^ and ft^^), Vdrn and Vigm in each cell. 

The time step is chosen by the minimum value among the following three time scales. The 
first is from the Courant condition given by 

St < n /^^^t^^", , , (18) 

maxd^j;! + Cs, \vy\ + Cs, \vz\ + Cg) 

where Ax is the cell size, Cg is the local sound speed, v^, Vy and Vz are the local fluid velocities 
and cfl is the Courant number, which we take as cfl = 0.6. The second time scale is imposed by 
cosmic expansion which requires that Aa/a < 0.02 within a single time step. The last time scale 
comes from the requirement that a particle moves not more than a fixed fraction of the cell size. 

The effect of the numerical resolution of these samples has been tested to higher order statis- 
tics(Feng, Pando &: Fang, 2003). It was found that two samples with different resolutions give 
about the same statistical result in the intermittency out to order eight. 

For statistical studies, we randomly sample 500 one-dimensional (ID) fields from the simulation 
results at redshifts z=4, 3, 2, 1, 0.5 and 0. Each ID sample, covering size L=25 h~-'^Mpc, contains 
192 data points, each one containing the mass density and peculiar velocity of dark matter, and 
the mass density and peculiar velocity of the IGM. 



1000 
100 
10 

1 

0.1 
0.01 

■S? 1000 
Q. 

100 



lz=0 



10 
1 

0.1b 
0.01 = 



T] I I MUM I I I INN I I MUM I I MIMI I I MUM I I I INN 




z=2 




I I MMMj I I MIIMj I I I IMMj I I MIMIj I I MIMIj I I MUM 

z=l 




M z=4 




al- I I I I I 

0.01 0.1 1 10 100 1000 



I I Ml I I I 

0.01 0.1 1 10 100 1000 



Fig. 1. — he relation between Pigm{x) and pdm{x) for redshifts z = 0, 1, 2 and 4. The data consists 
of 2± 7000 randomly drawn points from the simulation sample 

As first evidence of the discrepancy between the IGM and dark matter, Fig. 1 gives the relation 
between and Pigm at redshifts 4, 2, 0.5 and 0. If p^m traces pigm point-by-point, we should 
have a tight correlation along the line pdni = Pigm- However, Fig. 1 shows that although the 
relation p^^ = pigm is correct on average, the data points are significantly scattered around the 
line pdm = Pigm in the entire density range from p^m ~ 0.01 to 10^. 

The scatter of the IGM vs. dark matter has been noted in SPH simulations (Gnedin Sz 
Hui 1998). They found that the IGM density contrast is significantly different from that of dark 
matter. They correctly pointed out that a pure N-body simulation would fail to reproduce the 
IGM distribution without mimicking the dynamical effect of the gas. They explain the scatter as a 
dynamical effect of gas pressure. The pressure of the IGM with temperature ~ 10^ K is higher at 
higher redshifts, as the mean density of the gas is higher at earlier times. However, Fig. 1 shows 
that the scatter is smaller at higher redshifts. This would seem to contradict the gas pressure as an 



explanation for the scatter. On the other hand, the higher the nonhnearity, the higher the Reynolds 
number [eq.(6)], and the higher the discrepancy. Therefore, the gravitational nonlinear evolution 
of the IGM and dark matter system provides a plausible explanation of the scatter of Fig. 1. 



4. One-point Statistics 



4.1. One-point variables w^ith the DWT decomposition 



To calculate the one point distribution, we use the discrete wavelet transform (DWT). The 
scaling functions of the DWT analysis serve as the sampling window function. For the details of 
the mathematical properties of the DWT see Mallat (1989a,b); Meyer (1992); Daubechies (1992), 
and for cosmological applications see Fang & Thews (1998), Fang & Feng (2000). 

Let us briefly introduce the DWT-decomposition for a random field. Consider a 1-D density 
fluctuation d{x) on a spatial range from x = to L. We divide the space into 2^ segments labeled 
hy I = 0, 1, ...2-' — 1 each of size . The index j is a positive integer and gives the length scale 
The larger the j is, the smaller the length scale. Any reference to a property as a function 
of scale j below must be interpreted as the property at length scale L/2^ . The index / represents 
position and it corresponds to the spatial range lL/2^ < x < {I + 1)L/2K Hence, the space L is 
decomposed into cells 

The discrete wavelet is constructed such that each cell {j, I) supports a compact function, the 
scaling function <pj^i{x), which satisfies the orthonormal relation 



where S is Kronecker delta function. The scaling function (f)j^i {x) is a window function on scale j 
centered around the segment I. 

For a field F{x), its mean in cell {j, I) can be estimated by 




(19) 




So F{x)(j)j^i{x)dx 
Jo 4'j,ii^)dx 



Jq (l)j^i{x)dx 



1 




(20) 



where ef, 



is called scaling function coefficient (SFC), given by 




(21) 



Thus, a 1-D field F{x) can be decomposed into 



2^-1 




(22) 
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The term 0{> j) in eq.(22) contains only the fluctuations of the field F{x) on scales equal to and 
less than . This term does not have any contribution to the window sampling on scale j. Thus, 
for a given j, the one-point variables Fj^i or ej'j (/ = 0, 1...2-' — 1) give a complete description of the 
field F{x) smoothed on scale L/2K As one-point variables the ej^ are similar to the measure given 
by count-in-cell technique. However, the orthonormality eq.(19) insures that the set of Fj^i or t'^\ 
do not cause false correlations. In the calculations below, we use the Daubechies 4 (D4) wavelet 
(Daubechies, 1992). 
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Fig. 2. — The relation between pj^i-igm and Pj^i-dm on scales j = 5, 6, 7 and 8, at redshifts z = 0. 
The comoving scale corresponding to j is 33/2-^ h Mpc~^. The data consists of ~ 7000 randomly 
drawn points from the simulation sample 

Wc first calculate the one-point distributions of variables Pj^i-igm and Pj^i-dmi which are the 
one-point variables given by eq.(20) replacing F{x) by the density fields Pigm{x) and Pdmi^)- We 
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Fig. 3. — One-point distributions of the IGM density Pji-igm/ {p^ji^igm)^^'^ at redshift z = 1 and on 
scales j = 5, 6, 7, 8. The comoving scale corresponding to j is 33/2-' h~^ Mpc. 



then plot in Fig. 2 pj^i^igm vs. Pj^i-dm for j = 5, 6, 7 and 8, corresponding to comoving length 
scales 1.03, 0.516, 0.258 and 0.129 h~^ Mpc, respectively. Fig. 2 shows that the scatter around the 
line Pj^i-igm = Pj,i-dm is little smaller for smaller j. That is, the discrepancy between the IGM and 
dark matter is smaller on larger scales. Nevertheless, the discrepancy is still substantial on scale 
j = 5 or 1.03 h~^ Mpc, which is larger than the Jeans length of the IGM. 

When the "fair sample hypothesis" (Peebles 1980) holds, each set of 2^ {I = 0...2-'~^) points 
form an ensemble from which the one-point distribution can be studied. Figure 3 gives the one-point 
distribution of Pj^i-igm on scale j = 5, 6, 7 and 8 at z = 1. These distributions are generally non- 
Gaussian, having a longer tail at higher j, i.e., smaller scales. These distributions show significant 
j-dependence. The distribution at j = 8 has a tail longer than 9 times the variance, while the 
tail of the distribution of j = 5 is only about four times of the variance. We also see from Fig. 
3 that the distribution has large change from j = 5 to j = 6, and from j = 6 to j = 7, but a 
smaller change from j = 7 to j = 8. This is because the scale j = 8 is 0.129 h~^ Mpc which is 
close to the Jeans length. On these scales the perturbations in the IGM field are weak. We show 
the redshift-evolution of the pj^i-igm one-point distributions in Fig. 4. The one-point distributions 
also undergoes significant redshift evolution. However, a long tail PDF is already pronounced at 



] 




Fig. 4. — The redshift-dependence of the one-point distribution of the IGM density 
Pji-igm/ {Pji-igm)^^'^ ™ scales j = 6 and 8 corresponding to comoving length scales 0.52 and 0.13 
h"i Mpc 



5. Statistical discrepancy between the IGM and dark matter 

5.1. One-point distributions of velocity fields 

To demonstrate the statistical discrepancy addressed in §2, we first analyze the one-point 
distributions of the 1-D velocity fields Vigmi^) and Vdmi^)- For the linear solution vigmi^) = Vdm{^) 
[eq.(l)], we have Vj^i^igm = Vj^i^dm-, which are the one-point variables given by eq.(20) replacing 
F{x) by the velocity distribution. It is clear from Figure 5 that the one-point distributions of 
Vj^i-igm and Vj^i-dm at redshift z = are very different on all scales, j = 5, 6, 7 and 8. This 
difference is smaller at higher redshift. At z = 4, the one-point distributions of both Vj^i^igm and 
Vj^i-dm are basically the same on all scales. 

The nature of the statistical difference between the vj^-igm and vj^i-^dm one-point distributions 
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Fig. 5. — One-point distributions of the velocity fields of the IGM vji-igm and dark matter Vji-igm 
on scales j = 5, 6, 7, 8, and for redshifts z = (5a), 1 (5b), 2 (5c), 4 (5d). 



is not only quantitative, but also qualitative. The Vdm{^) one-point distributions are small devi- 
ations from a Gaussian PDF on all scales and all redshifts considered. This result is consistent 
with previous studies using N-body simulation samples (e.g. Yang et al 2001). On the other hand, 
the one-point distributions of the IGM velocity field shown in Fig. 5 are generally exponential at 
redshifts z < 2. Even when the scale is as large as j = 5 or 1 h^"*^ Mpc, the distribution of Vigm{'^) 
is still exponential (Fig. 5a). This emphatically shows that the statistical discrepancy between the 
IGM and dark matter developed with dynamical evolution. 

The dynamical equation (3) for the dark matter looks very similar to eq.(4) for the IGM. So 
why do the two one-point distributions have such different shapes? The reason is that the cj) term 
of eq.(4) is truly an external force, as the gravity potential is independent of the IGM density and 
velocity, dependent only on the dark matter. Therefore, when its Reynolds number is large at lower 
redshifts, shocks or Burgers' turbulence will develop in the IGM field due to the external driving 
force and the IGM velocity field will be highly non-Gaussian. On the other hand, the dark matter 
mass density in eq.(3) is not independent of the gravitational potential and there is no external 
driving force. Thus, when the gravity potential ^ is Gaussian in weakly non-linear evolution, the 
dark matter velocity PDF will approximately be Gaussian too. 



-16- 




Fig. 6. — One-point distributions of the difference Vji-dm — Vji-igm on scales j = 5, 6, 7, 8, and for 
redshifts z = (6a), 1 (6b), 2 (6c), 4 (6d). 

The statistical discrepancy of the velocity fields can also be seen with Fig. 6 which gives 
the one-point distribution of the difference Av^y = Vj^i-dm — ''^j,l-igm on scales j = 5, 6, 7 and 
8, and redshifts 0, 1,2 and 4. These one-point distributions are non-Gaussian on all scales and 
redshift considered. Although the one-point distributions of Vdmi^) and Vig„i{x.) at redshift z = 4 
are not very different [Fig. 5d], their differences [Fig. 6d] are highly non-Gaussian. Therefore, 
the discrepancy between vj^i-dm a-^d Vj^i^igm is not due to noise, but arises from the non-linear 
evolution of the IGM fluid. 



5.2. One-point distributions of density fields 

Figure 7 presents the one-point distributions for pj^i^igm and pj^i-dm on scale j = 7 and redshifts 
z = 1, 2 and 3. The horizontal axes are the variance-normalized densities, Pj,i-igm/ {p] i-igm)^^'^ 
and Pj,i-igm/ {Pj i-dmy^^ ■ Unlike the velocity fields, both the IGM and dark matter show highly 
non-Gaussian behavior even at redshift z = 3. Additionally, the two curves in Fig. 7 are not 
significantly different. That is, the mass density one-point distribution is not as sensitive to the 
statistical discrepancy as the velocity one-point distribution. This is probably because the statistical 




Fig. 7. — One-point distributions of the IGM density Pji-igm/ {Pji-igm)^^'^ and dark matter 
Pji-dm/ {p'ji^d-m)^^'^ redshift z = 1, 2 and 3 and on scales j = 7, or comoving length scale 
0.26 h-i Mpc. 

discrepancy between a "passive substance" and underlying field is significant only for non-conserved 
quantities (temperature, velocity etc.) (Shraiman & Siggia 2001.) Nevertheless, one can see from 
Fig. 7 that the IGM distributions are always a little higher than the corresponding distribution for 
the dark matter at the range around Pj,i-igm/ {p'j i^igm)^^'^ — 1- This indicates that the tail of the 
IGM one-point distribution should be shorter than that of dark matter. Therefore, we can expect 
that the high order moments of the one-point distributions will more clearly show the discrepancy 
between the IGM and dark matter. 



In this paper, we use only the second order moment. The correlation between pj , 



-igm 



and 

Pj,i-dm can be measured with the ratio Idm, ligm [eqs.(lO), (11)], or / [eq.(12)]. Using the DWT 
one-point variables Pjj—igm 

and Pj^i-dm, we can redefine, respectively, the ratios Idm and ligm as 
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Fig. 8. — The ratio Ij at redshifts z = 0, top left panel, to redhsifts 0.5, 1, 2, 3, 4, (moving left to 
right) and on scales j = 5, 6, 7, 8, which corresponds to the comoving scales 33/2^ h'^ Mpc. 



Prom the one-point distributions of Fig. 7, we have {p\i_dm) — {Pj,i-igm)- Thus, if Pj^i-igm per- 
fectly traces Pj,i-dm, wc have correlation {pj,i-igmPj,i-dm) ^ {p],i-igm) - {P%i-dm)^ therefore, 
Ij-dm - Ij-igm - 0. If Pj,i-igm is fully independent of Pj,i-dm, we have Ij-dm - Ij,-igm - V^. 

Figure 8 plots Ij-dm and Ij-igm for z = 0, 0,5, 1, 2, 3 and 4 on scales j = 5, 6, 7 and 8. The 
error bars are the 67% confidence level of the 500 1-D samples. Figure 8 shows that Ij is always 
larger than 0, equal to ~ 0.1 — 0.2 at redshift > 2, and ~ 0.3 — 0.6 at redshift < 2, and weakly 
dependent on scales. That is, the cross correlation at redshift < 2 is 

{Pj,l-igmPj,l-dm) {Pj,l-igmPj,l-dm) ^ r, 7n fic; fOr.\ 

\Pj,l-igm/ \Pjd-dml 

Therefore, the IGM is not a perfect tracer at low redshift and on scales larger than the Jeans length. 
These results indicate that 15 -30% of the two density fields are off-correlated. This is consistent 
with the analysis in §2.2. 
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Fig. 9. — One-point distributions of the difference Aj^i/ {A'j i)^/"^ > for redshift z = 1, 2 and 3, 
and on scale j = 7, or 0.26 Mpc. 

5.3. One-point distribution of A(5(x) 

Tlie DWT one-point variable for the density difference AS(x) [eq.(2)] is given by 

= J A5{x)(t>j,i{x)dx. (26) 

Since (t)j,i{x) is localized in cell the integral J (pj^i(pj^i/(pj^i"dx is generally non-zero only for I = 
I' = I" (for Haar wavelet, this is exactly true. For other wavelets, the integral with I / I" is much 
less than that oil = 1' = I"). Thus, for two fields A{x) and B{x), we have J" A{x)B{x)il)ji{x)dx 2± 
dj^fl^fli where e^^ and are the WFCs of A{x) and B{x). The factor gj is 



9j = 

which is the factor gw used in eq.(16). 

With Aj^i, eq. (15) can be rewritten as 



J il^%{x)dx. 



(27) 



(28) 
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where 

r,j = I dx'^ji{x'){V . Av{x)). (29) 

and 

Xj = gj I x'^ji {x') (V -Ydm). (30) 
As discussed in §2.3, the PDF of Aj ; should be long tailed. 

The normalized one-point distributions of Aj^i/ {A'j i)^^'^ for j = 7 and redshifts 1,2, and 3 
arc plotted in Fig. 9. All the tails of the distributions are remarkably long. They have 6 a 
events at 2; < 3 with probability ~ 0.005, while for a Gaussian field, the probability of 6-cr is 
< 10~^. This shows again that the statistical discrepancy between the IGM and dark matter is 
not due to processes like Jeans diffusion or Gaussian noise. Fitting the tails with a power law 
(Aj^;/(A|^)^/^)~^, we have (3 <1. This result is consistent with the dynamical eq.(28). 



6. Conclusion 

Using the dynamical equations governing the IGM, we have shown that the decoupling of 
important statistical properties of the gas component from the underlying dark matter field on 
scales larger than the Jeans length is inevitable. The one-point distribution of the IGM peculiar 
velocity field is found to be substantially different from that of dark matter at redshift < 2. Although 
the one-point distribution of the IGM mass density field shows only a little deviation from that of 
the dark matter, the density fields of the IGM and dark matter are not point-by-point correlated, 
and a significant part of the two fields is off-correlated on all redshifts z < 4. 

It is difficult to explain the discrepancy as an effect of Jeans smoothing, as the discrepancy is 
still evident even when the fields are smoothed on scales > 1 h~^ Mpc (j = 5). This is much larger 
than the Jeans length of the IGM with temperature ~ 10^ K at the mean density. The Jeans length 
of the IGM may be > 1 h~^ Mpc in an area with high temperatures ~ 10^ K, but not one with 
high density. Yet, even in this case Jeans smoothing does not adequately explain the statistical 
discrepancy, because the difference between the density distributions of the IGM and dark matter 
is highly non-Gaussian. 

The decoupling of the IGM from the dark matter mass field is important for problems in 

large scale structure formation. Some observations have already implied that the IGM and dark 
matter have decoupled. For instance. X-ray measurements find no evidence for the baryon fraction 
of clusters to be equal to the universal value from cosmological nucleosynthesis (White & Fabian 
1995; David 1997; Ettori, Fabian & White 1997; White, Jones & Forman 1997; Ettori & Fabian, 
1999). This result violates the similarity the IGM and dark matter if typical galaxy clusters formed 
from linear fields on scales of a few ten (comoving) Mpc, which is much larger than the Jeans length 
of the IGM. Moreover, the X-ray observed luminosity-temperature relation for groups and clusters 
is found to be inconsistent with the prediction given by the similarity between the IGM and dark 
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matter. The high entropy floor observed in nearby groups and low mass clusters directly violates 
the dynamical self-similar scaling (Ponman, Cannon, & Navarro 1999; Lloyd-Davies, Ponman, & 
Cannon 2000). Although, at present we cannot attribute all these discrepancies to the decoupling 
shown in this paper, we now understand that the dynamics of the IGM and dark matter fields 
will lead to a qualitatively different evolution of those fields. While we did not consider stellar 
formation and its feedback to the IGM in our present simulations, we believe that the statistical 
discrepancy will still be a common feature even when these effects are considered because the reason 
for discrepancy we have uncovered arises from the nonlinear evolution of the random fields of the 
IGM and dark matter. 

Finally, we should mention the effect of the size of the simulated sample. Since the simulation 
box is only 25 h~^ Mpc, the use of the fair sample hypothesis may not be appropriate when 
considering perturbations on scales larger than this. We can estimate the effect of these larger 
perturbations by adding long- wavelength modes (Tormen &; Bertscinger 1996). Since the IGM 
and dark matter fields are linear or quasi-linear on scales larger than 25 h~^ Mpc, the effect of 
long-wavelength perturbations can be analyzed by adding a displacement to each particle in the 
simulation box. The displacement is given by the linear field or Zeldovich approximation of the 
field consisting of modes on scales larger than 25 h^-*^ Mpc. As was discussed in §1, there is no 
discrepancy between the IGM and dark matter in the linear or quasi-linear regime. Thus, the 
displacement for the IGM is the same as that of the dark matter meaning that long wavelength 
modes are not a source of the discrepancy between the IGM and dark matter. 

We thank Drs. C.W Shu and P. He for their help during the preparation of this revised version. 
Thanks also to our referee for his (her) comments. LLF acknowledges support from the National 
Science Foundation of China (NSFC) and National Key Basic Research Science Foundation. 



A. Hydrodynamic equations for dark matter fields 

Let us consider a flat universe having cosmic scale factor a{t) oc t^/^ and dominated by dark 
matter. In a hydrodynamic description, the dark matter is described by a mass density held 
Pdmip^it) and a peculiar velocity held Vdm(x, t), where x is the comoving coordinate. The field is 
described by the equations of continuity, momentum and gravitational potential as (Wasserman 

= (Al) 

(A2) 
(A3) 

The mean density is pdm{t) = l/Gi^Gt^ oc a~^. The gravitational potential ^ is zero (or constant) 
when the density perturbation 6dm = 0- The operator V acts on the comoving coordinate x. 



+ -V ■{! + 5dm)'Vdm 

at a 

—Q^ 1- (Vdm • V)Vd^ = - 



- 22 - 



For growth modes in the perturbations, velocity is irrotational. We can then define a velocity 
potential by 

Vdm = "^^Pdm- (A4) 

a 

The momentum equation (A2) can then be rewritten as 
This is the Bernoulli equation. 



B. Hydrodynamic equations for the IGM 

As usual, the IGM is assumed to be ideal fluid with polytropic index 7 = 5/3. The hydrody- 
namic equations of the IGM are (Peebles 1980) 

^ + iv . (1 + Sigm)-Vigm = (Bl) 
+ i^igm ■ V)^^^gm = ^Vp - V<A (B2) 

1 1 1 

— + 5-S + -V • (Svigm) = V • (pViyrn) PigmVigm ' V(/> - A^ad, (B3) 

01 a a a a 

where Pigm, ^igm, £ and p are, respectively, the mass density, peculiar velocity, energy density 
and pressure of the IGM. The term Aj-^d in Eq.(2) is given by the radiative heating-cooling of the 
baryonic gas per unit volume. The gravitational potential in eqs.(B2) and (B3) can still be given 
by eq. (A3). That is, the gravity of the IGM is negligible. The evolution of the IGM mass field 
Pigm(x, t) is governed by the gravity of dark matter only. 

The hydrodynamic equations for the IGM, eqs.(Bl)-(B3) can be written in the form of con- 
servation laws for mass, momentum, and energy in a comoving volume as 

— df^ ~^ a "'^P^^^^ ^ -aa^pvi - a'^pV(f> (B5) 



''J 

do?£ 1 d 



dt ^ ad^-^^^^^'^'"'^ " -2aa^f - a^pv • V<^ - o'Kad- (B6) 
In eqs.(B4)-(B6), we dropped the subscript igm for simplicity. 

To sketch the gravitational clustering of the IGM, it is not necessary to consider the details of 
heating and cooling. Thermal processes are generally highly localized, and therefore, it is reasonable 
to describe all thermal processes by a polytropic relation p(x, t) oc p]g^{'x., t). Thus eq.(B2) becomes 

— + (v.,„.V)v.^ = - — ^^-^-V* (B7) 
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where the parameter ji is the mean molecular weight of the IGM particles, and rup the proton 
mass. Here, we don't need the energy equation and the IGM temperature evolves as T oc p'^~^, or 
T = TQ(l + 6igmV-\ 

Eq.(B7) differs from eq.(A2) only by the temperature-dependent term. If we treat this term 
in the linear approximation, we have 

- ^ {'^V^an^? " ^'^\^gm = <P, (B8) 

where ipigm is the velocity potential for the IGM field defined by 

Vigm = -^Vipigm- (B9) 

The coefficient u is given by 

fimp{dlnD{t)/dt)' ^ ' 

where D(t) describes the linear growth behavior. The term with v in cq.(B8) acts like a viscosity 
(due to thermal diffusion) characterized by the Jeans length kj = {a^ /t^){ump/-fkBTQ). 

The unperturbed solutions of the density and velocity fields of both dark and baryonic matter 
are = {Q.i,/^dm)p ^ a-^d v^ = v = 0, where and ^dm a-i's, respectively, the density 
parameters of the IGM and dark matter. Therefore, the linearlization of eqs. (Bl) and (B7) yields 

^ + Iv • v,,„ = (Bll) 

where the mean temperature T oc 'pb^~^ oc a'~^^'^~^\ In Fourier space, we have 

a — di -^^^i9^y^^'^) = "^T^GpSdmiKt) (B13) 



^}^1EI^ + \g^O,,t) = --l^<5,^^(k,t) + ^^6dm{Kt), (B14) 



where Vj^mCk,*) = i^'^igmO^jt), and the Jeans wavenumber kj = {a^ /t'^){fimp/jkbT) oc a^'' ^. If 
7 = 4/3, kj is time-independent. 

In solving equations (B13) and (B14), we consider only the growth mode of the perturbation 
of dark matter, i.e. 5(k,t) oc a. In the case of 7 = 4/3, the solution of eqs.(B13) and (B14) is (Bi, 
Borner, Chu 1993) 
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where e = (1 — 4/c^/9A;j)^/^, and constants ci and C2 depend on the initial condition 5igrre(k, 0) and 
VigmO^,0)- Therefore, regardless the initial condition of IGM, after a long evolution we have 

SigmiKt) = SdmiKt), VigmiKt) = Vd„(k,t), if < fcj. (B16) 

These solutions mean that the initial conditions of the IGM are unimportant. The IGM will, in 
the end, follow the same trajectory as the dark matter on scales larger than the Jeans length. 

The solutions of linear equations (Bll) and (B12) have also been found by using assumptions 
of the IGM thermal processes other than that used in eqs.(B13) and (B14) (e.g. Nusser 2000, 
Matarrese Sz Mohayee 2002). A common feature of these solutions is 

Sb(k, t) = {1 + decaying terms)(5(k, t) + decaying terms, if fc <C kj. (B17) 

The initial conditions of the IGM field affects only the decaying terms, and therefore, the linear 
solutions eq.(B16) hold in general regardless specific assumptions of IGM thermal processes. 



C. PDF of Aj^i in eq.(27) 
Using dr = dt/a, eq. (28) can be rewritten as 

^ = -x,A,,+n,. (ci) 

If the stochastic forces Xj and rjj are Gaussian, then 

(A,) = A,- (C2) 
([A,(r)-A,][A,(rO-A,]) = 2D^jS{t - r') 

iVj) = 
{r]j{r)r]j{T')) = = 2D,,,-5(t - r') 

The Fokker-Planck equation corresponding to eq.(Cl) is (Venkataramani et al. 1996) 

where P(Aj ;,r) is the PDF of A^ ;, and the flux j(Aj ;,r) is given by 

j{^j,i,r) = (-A, + D^j) Aj,i P{Aj,i,T) - ^[Dx,,\ Aj,i p + I),,,]P(A,y, r). (C4) 

For stationary solution, dP{Aj^i,T)/dT = 0, and therefore, dj(Aj^i,T)/dAj^i = 0. We have 
j(Aj ;) = const. However, when A^y is very large, we should have i(Aj ;) = 0, and therefore 
j{S^) = 0. Thus, we have equation as 

+ Dxj) Aj,i P{Aj,i, r) - ^[DxJ Aj^i 1^ + D^j]P{Aj,i,t) = (C5) 
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The solution of eq.(C5) is 



(C6) 



where C is a normalization constant. Th PDF is then 




(C7) 



where 




s = 



(C8) 



and 




+ 1. 



(C9) 
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